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We theoretically study nonequilibrium spin transport in a superconducting wire connected 
by tunnel junctions to two ferromagnetic metal wires, each of which serves as an injector or de- 
tector of spin-polarized electron current. We present a set of Boltzmann equations to determine 
nonequilibrium quasiparticle distributions in this system, and obtain an analytical expression 
for the nonlocal spin signal in the case of small injection current. It is shown that the quasipar- 
ticle distribution in the ferromagnetic metal for detection strongly affects the magnitude of the 
spin signal. At low temperatures, since nonequilibrium quasiparticles created by the tunneling 
from the superconductor dominate thermally excited ones, the spin signal becomes independent 
of temperature. This explains the convergence of the spin signal with decreasing temperature 
observed in a recent experiment by Poll et al. 
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1. Introduction 

Experimental studies on spin injection and detection 
in a normal metal have attracted considerable atten- 
tion recently in the field of spintronics.^"^ More than 
two decades ago, Johnson and Silsbee^ performed the 
first experiment on this subject by using a large normal 
metal sample with two electrodes made of a ferromag- 
netic metal, where each electrode serves as a spin injector 
or detector. Spin-polarized electrons created near the in- 
jector diffuse in the normal metal, and spin imbalance is 
transmitted to the detector if spin-fiip scattering does not 
suppress it. They found an evidence of spin imbalance by 
measuring an open-circuit voltage induced at the detec- 
tor. Several experiments using devices in the mesoscopic 
regime have been reported to date.^^^ The most popular 
device consists of a thin normal metal wire connected to 
a few ferromagnetic metal wires. In this system, we sup- 
ply injection current /jnj with spin polarization Pgpin into 
the normal metal from one of ferromagnetic metals and 
measure an open-circuit voltage between another ferro- 
magnetic metal and the normal metal. Let Vp (T4p) be 
the open-circuit voltage when the magnetizations of the 
two ferromagnetic metals are parallel (antiparallel). We 
are interested in the nonlocal spin signal defined by 



ap 



(1) 



which crucially depends on the spin diffusion length Agf 
and the distance d between the injection and detection 
points. In the case where the normal metal and the two 
ferromagnetic metals are connected by tunnel junctions, 
the spin signal is given by"^' ^ 



Pspi„i?Ne 



(2) 



where i?N = pNAsf/j4N with pN and being the resis- 
tivity and the cross-sectional area of the normal metal, 
respectively. 

Spin injection and detection in a superconductor is 
also attracted considerable attention. ^"^^ Our primary 



interest focuses on how the spin signal is modified by the 
transition to the superconducting state. Takahashi and 
Maekawa^ studied this problem and predicted that 

^ ,i?Ne"^, (3) 



-^spin 



2/o(A)^'p-^ 



where /o(A) = l/(exp(A/T) + l) with the superconduct- 
ing energy gap A and temperature T . This indicates that 
Pspin exponentially increases with decreasing T. They 
claimed that this modification is caused by the increase 
of spin resistivity due to the opening of the energy gap 
A. The increase of i?spin with decreasing T has been 
successfully observed in the recent experiment by Poll 
et al}^ However, there remain a few points to be clari- 
fied. We focus on the following two points. Firstly, Taka- 
hashi and Maekawa implicitly assume in their derivation 
of eq. (3) that spin imbalance in a superconductor can 
be described by a shift of spin-dependent chemical poten- 
tial. This assumption cannot be justified at low tempera- 
tures, where energy relaxation due to phonon scattering 
is not strong. Secondly, Poli et al. observed convergence 
of i?spin with decreasing temperature. This behavior can- 
not be explained by eq. (3). 

In this paper, we theoretically study nonequilibrium 
spin transport in a hybrid system consisting of a super- 
conducting wire and two ferromagnetic metal wires. Each 
ferromagnetic metal is connected by a tunnel junction 
to the superconductor, and serves as an injector or de- 
tector of spin-polarized quasiparticles. We present a set 
of Boltzmann equations governing nonequilibrium quasi- 
particles in this system. We focus on the case of small 
injection current at low temperatures, and obtain not 
only the quasiparticle distribution in the superconduct- 
ing wire but also that in the ferromagnetic metal wire 
for detection. On the basis of the resulting nonequilib- 
rium distributions, we derive an analytical expression for 
the nonlocal spin signal. It is shown that although the 
spin signal originates from spin imbalance transmitted to 
the detection junction, its magnitude is not solely deter- 
mined by the spin imbalance but is strongly affected by 
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the quasiparticle distribution in the ferromagnetic metal. 
We observe that when T is higher than a crossover tem- 
perature Tcross, the spin signal exponentially increases 
with decreasing T reflecting the reduction of thermally 
excited quasiparticles in the ferromagnetic metal. At low 
temperatures below Tdoss. however, the magnitude of the 
spin signal is determined by noncqniUbrium quasiparti- 
cles created by the tunneling from the superconductor 
instead of thermally excited ones, and the spin signal be- 
comes independent of T . This explains the convergence 
of the spin signal with decreasing T observed by Poll et 

In the next section, we present a set of Boltzmann 
equations to describe nonequilibrium quasiparticle dis- 
tributions in the hybrid system consisting of a supercon- 
ducting wire and two ferromagnetic metal wires. In §3, 
we obtain nonequilibrium quasiparticle distributions in 
this system by solving the set of Boltzmann equations, 
and derive an analytical expression of the spin signal on 
the basis of the resulting quasiparticle distributions. In 
§4, we compare our theoretical result with the recent ex- 
perimental result. We set ?l = fee = 1 throughout this 
paper. 

2. Formulation 

Let us consider the hybrid system consisting of a su- 
perconducting wire and two ferromagnetic metal wires 
(see Fig. 1). We assume that the superconductor is con- 
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Fig. 1. (a) Schematic picture of the model system consisting of a 
superconducting wire (S) and two ferromagnetic metal wires (F). 
The left (right) ferromagnetic metal serves as an injector (detec- 
tor) of spin-polarized electron current, (b) Spatial coordinates 
used in the text. 



nected by a tunnel junction to each ferromagnetic metal. 
The left and right junctions serve as spin injector and de- 
tector, respectively. We adopt a simple one-dimensional 
model for this device assuming that the superconductor 
and the ferromagnetic metals are very thin. We intro- 
duce the X axis in the superconductor on which the left 
and right junctions are located at a; = Xinj and x = Xdet, 
respectively, and the y axis in the left (right) ferromag- 
netic metal on which the injection (detection) junction 



is located a.t y — yinj (j/dot)- We denote by d the separa- 
tion between the two junctions. That is, d = Xdet — Xiny 
We inject spin-polarized current into the superconduc- 
tor by applying a bias voltage l^nj across the injection 
junction, and measure an induced open-circuit voltage 
Vdet across the detection junction under the condition 
that net current flow vanishes between the superconduc- 
tor and the ferromagnetic metal for detection. We sim- 
ply assume that the spin polarization Pgpin of the injec- 
tion current is proportional to the difference between the 
density of states Nf-[ for up-spin electrons and Npi for 
down-spin electrons. The spin polarization is expressed 



as 



FT 



Fi 



■ spin 



(4) 

iVpt + Nfi ^ ' 

We assume that spin relaxation in the superconductor is 
caused by spin-flip scattering due to spin-orbit interac- 
tion as well as magnetic impurities. 

To present an expression for the tunneling current 
across each junction, we consider nonequilibrium quasi- 
particle distributions in the superconductor and the fer- 
romagnetic metals. We first introduce the quasiparticle 
distribution function S'flct in the left ferromagnetic metal 
for injection, where a =Tji is the spin variable. We as- 
sume 5fl<t(2/) e) = /o(e — e^nj) with the Fermi-Dirac dis- 
tribution function /o(e). Here and hereafter, we measure 
quasiparticle energy from the chemical potential of the 
superconductor not only in the superconductor but also 
in the ferromagnetic metals. We next introduce the quasi- 
particle distribution fimction ggcr in the superconductor. 
In terms of four distribution functions /l-, /t+ 

and /t- for nonequilibrium quasiparticles, we express it 
asi2-i6 

5'ST(2;,e) = /o(e) + /L-F(a;,e) +/L-(x,e) 

+ /T+(.x,e)+/T-(x,e), (5) 

584(2;, e) = /o(e) - /L-K(x,e) +/L_(x,e) 

+ /T+(x,e)-/T-(x,e). (6) 

The four distribution functions satisfy 

.fL,T+(.T, -e) = /L,T+(a;, e), (7) 

/L,T-(.x',-e) = -/L,T-(x,e). (8) 

Note that /l+ describes spin imbalance, while /t+ de- 
scribes charge imbalance. The other two functions 
/l_ and /t- describe total energy imbalance and en- 
ergy imbalance between up-spin and down-spin quasipar- 
ticles, respectively. Finally, we introduce the distribution 
function g^Bjj in the right ferromagnetic metal in which 
nonequilibrium quasiparticles appear due to quasiparti- 
cle tunneling from the superconductor. We express it as 



9FR^{y, e) = /o(e - eVdet) + fpaiv, e)- 



(9) 



We hereafter assume that the magnitude of the energy 
gap A is unaffected by spin injection everywhere in the 
superconductor. This allows us to consider fh±{x, e) and 
fT±{x,e) only for |e| > A. The nonequilibrium distribu- 
tion functions /l±, /t± and /pcr are governed by Boltz- 
mann equations which we present below. 
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The tunneling current at the injection junction is given 



by 



-Jl(Mnj,T), 



(10) 



where Rinj is the tunnel resistance of the injection junc- 
tion and 

1 f°° 

Ji {V, T) = ^j deTVi (e) (/o (e - eV) - /o (e + e^)) 

(11) 

with Ni being the normalized density of states in the su- 
perconductor, given by -/Vi(e) = |e|/-\/e^ — for |e| > A 
in the BCS limit. In deriving eq. (10), we have ignored 
small contributions arising from noncquilibrium quasi- 
particles in the superconductor. The tunneling current 
between the superconductor and the right ferromagnetic 
metal for detection is expressed as 

/dct(Kict) = /q(Kict) + /F(V;ict) - /s(Kict), (12) 

where Jq is the ordinary tunneling current arising from 
thermally excited quasiparticles, while /p represents the 
contribution from noncquilibrium quasiparticles induced 
in the ferromagnetic metal. The third term Jg represents 
the contribution from spin and charge imbalances. They 
are expressed as 



/q(Klet) 



-Jl(Klet,T), 



/F(Kiet) = -^ rdeN,{e) 
e^tdet Jo 



(13) 



1 + K 



spin 



(/FT(2/det,e) + /FT(2/det, -e)) 



+ 



spin 



f 

Jo 



ifFiiVdet, e) + fpiiVdet, -e)) 



deiVi(e) 



(14) 



eRdet 

X (fspm/L+(a;det,e) + /T+(a;det,e)), (15) 

where i?dot is the tunnel resistance of the detection junc- 
tion. In eq. (15), the first term with /l+ represents the 
contribution from spin imbalance and is the origin of the 
spin signal, while the second term with /x+ represents 
that from charge imbalance. In deriving eqs. (14) and 
(15), we have assumed the parallel alignment of magne- 
tizations. The corresponding expressions for the antipar- 
allel alignment is obtained by reversing the sign of Pspin- 
To present Boltzmann equations for /l± and /t±) we 
introduce the Usadcl cquation^^ for the quasiclassical re- 
tarded Green's functions and 

ie/«(e) + A5^(e) - ^5«(e)/«(e) = 0, (16) 

Tm 

where Tm represents the magnetic impurity scattering 
time and we have assumed that the superconductor is 
spatially homogeneous. The spectral functions Ni, N2, 
Ri and i?2 are defined as 

g^(e) =iVi(e)+ii?i(e), (17) 
/^(e)=7V2(e)+ii?2(e). (18) 
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In terms of the spectral functions, the Boltzmann equa- 
tions are expressed as"'^^"'^^ 

Ds {Nl{e)-Rl{e))dlh+{x,e) 

-^{Nl{e)-Rl{e))h+{x,e) 

' SO 

-^(^i'(e)+i?i(e))/L+(x,e) 

+ PL+(ar,e)=0, (19) 
Ds {Nl{e) - Rl{e)) dlh-{x, e) + Pl-{x, e) = 0, (20) 
Ds {N^{e)+Ni{e))d'jT+{x,e) 
1 



^conv (^, 
r2 



■/T+(x,e)-f PT+(x,e) = 0, 



(21) 



Ds {N!{e)+Ni{e))dlh-{x,e) 
4 

3Tso 

4 

3Tm 
1 



(7Vf(e) + iV|(e))/T-(x,e) 
(iV^(e)-iV|(e))/T-(^,e) 
fT-{x,e) + PT-ix,e) = 0, 



(22) 



PL-(a;,e) 



where Ds is the diffusion constant, Tgo and Tconv are 
the spin-orbit scattering time and the charge imbalance 
conversion time, respectively, and Pl± and Pt± arc the 
injection terms which represent quasiparticle tunneling 
between the superconductor and the left ferromagnetic 
metal. The injection terms are given as^^'^° 

5{x - Xir.i)Ni{e) 



Jospin (/o (e - e^inj) - /o (e + eV-^j)) 

- 2(/L+(a;inj, e) + Pspin/T+(a;inj, e))] , (23) 

_ 5{x- a;inj)^i(e) 
4e27VsAsPinj 

/o (e + elinj) + /o (e - e^nj) - 2/o (e) 

- 2(/L-(a:inj, e) + Pspi„/T-(xi„j, e))] , (24) 
_ 5{x - Xini)Ni{e) 

462iVsAsPinj 

/o(e-eyi„j)-/o(e + eT4nj) 

- 2(Pspin/L-F(a;inj,e) + /T-F(a;inj,e)) , (25) 

_ 5{x - .i;inj).Vi(.r.c) 
4e27VsAsPinj 

npin(/o (e + eMnj) + /o (e - e^nj) - 2/o (e)) 



PT+{x,e) 



PT-(x,e) 



- 2(Pspin/L-(a;inj, e) + /T-(a;inj, e)) 



(26) 



where A^s and As are the density of states at the Fermi 
level in the normal state and the cross-sectional area of 
the superconducting wire, respectively. We can ignore 



4 J. Phys. Soc. Jpn. 



Full Paper 



Y. Takane 



/l± (a^inj , e) and /t± (a^inj , e) in these injection terms when 
the injection current is small. It should be noted that in- 
elastic phonon scattering has been ignored in the Boltz- 
mann equations because its role is not relevant at low 
temperatures, in which we are interested. We have also 
ignored very small contributions to /l± and /t± arising 
from the coupling with the right ferromagnetic metal for 
detection. 

We turn to quasiparticle distributions in the right fer- 
romagnetic metal for detection. We note that in obtain- 
ing Ip, the spin-dependence of fpa is not important as 
long as the spin polarization is small. This indicates that 
we need not consider complicated spin-dependent dy- 
namics of nonequilibrium quasiparticles. We thus define 

/FT(y,e) + /Fi(y,e) 



fF+{y,e) 



(27) 



and approximate the expression of Ip as 

1 r°° 

/F(Kiet) = ^— / deiVi(e) 
eKdct Jo 

X (/F+(?Aict, e) + fr+iVdet, -e))- (28) 

We present an appropriate Boltzmann equation for /f+ • 
We ignore roles of spin-flip scattering since the spin- 
dependence is not important for our argument. How- 
ever, we miist consider the energy relaxation process due 
to phonon scattering. The reason for this is as follows. 
Since quasiparticles in the ferromagnetic metal are in- 
duced by the tunneling from the superconductor with 
the energy gap A, their excitation energy is of the or- 
der of A and no quasiparticle is directly created in the 
subgap region. Quasiparticles in such a nonequilibrium 
situation inevitably experience the energy relaxation. We 
thus assume that fF+ obeys 

DpdyfF+iy, e) 7 — ^—rr^fF+{y, e) + Pp+iv, e) = 0, 

(29) 

where Dp is the diffusion constant averaged over spin di- 
rections and Te is the energy relaxation time with e— eVdet 
being the quasiparticle energy measured from the chemi- 
cal potential of the ferromagnetic metal. The source term 
Pf+ describing quasiparticle tunneling from the super- 
conductor is given by 

o / X (5(2/ - ydet)-^Vi(e) / , , . . 

= 2e^7VpApflaet '^"^'^ 

- fF+iVdet, e) + /L-(a;det, e) + h+ixdet, e)), 

(30) 

where Np = {Npi +Npi)/2 and Ap is the cross-sectional 
area of the ferromagnetic metal. For the expression of Te, 
we adopt 



de'ap{e,e') 



X ^coth 



2T 



tanh 



2r 



with 



ap{e,€') 



ap 



sign(e' — e) X (e' — e)^ 



(31) 



(32) 



where ap characterizes the strength of electron-phonon 
coupling. For |e| T, we approximately obtain 

1 



re(e) 



3 ' ' ■ 



(33) 



3. Spin Signal 

In this section, we solve the Boltzmann equations and 

obtain the spin signal defined in eq. (1) by evaluating Vp 
and Vap. Note that Vp (Kp) is the open-circuit voltage 
induced across the detection junction when the magne- 
tizations of the injector and detector are in the paral- 
lel (antiparallel) alignment. We determine Vp and V^ap 
by the condition of /det = 0. We assume that the mag- 
nitude of Vp and Vap is much smaller than A/e. How- 
ever, we do not assume Vinj A/e. We focus on the 
case where the injection current is so small that injected 
quasiparticles are populated only near the gap edge (i.e., 
|e| ~ A). In this case, /t-i- and /t- quickly relaxes be- 
cause the conversion time becomes very short near the 
gap edge.^^' ^^'^-'^'^^ Therefore, we ignore /t+ and /t- 
in the following argument. Furthermore, the smallness of 
the injection current also allows us to ignore /L+(a;inj, e) 
and /L-(a;inj, e) in the injection terms given in eqs. (23) 
and (24). 

We first assume that magnetic impurities are absent 
(i.e., = 0) and define the spin-flip scattering time Tgf 

as 



1 

Tsf 



4 
3ts( 



(34) 



In this case, the spectral functions for |e| > A are simply 
given by 



R2{e) 



Ve^- A^' 
sign(e)A 



(35) 



(36) 



Ve^- A^' 

and A^2(e) = Ri{e) = 0. This indicates that N^{e) - 
i?|(e) = 1. We first obtain /L+(a;detje) by solving 
eq. (19). Note that /l+(x, e) decays exponentially as a 
function of \x — Xi^jl and this decay is characterized by 
the spin-diffusion length given by Agf = V-DsTgf. We ob- 
tain 

/L+(a^det,e)^ ,^,,^Yo'' ^ E+{e,V^^)e-^ (37) 



SemsAsRinsDs 



with 



E+(e, Vini) = N,{e) (/o (e - e^nj) - /o (e + eVmi)) • 

(38) 

Next, we obtain /L_(a;det,e) which is necessary to ob- 
tain /F+(ydot:e)- A special care must be paid in solv- 
ing eq. (20) since no relaxation process is included in 
this equation. The relaxation of /l- is mainly caused 

by the phonon-mediated recombination process, which is 
described by adding the following nonlinear term^° 

/L-(x,e) = -4 / de'as{e,e') 



X {m{e)m{e') - i?2(e)i?2(e'))/L-(a;,e)/L-(x,e') 

(39) 



J. Phys. Soc. Jpn. 



Full Paper 



Y. Takane 5 



to eq. (20). Here, (Ts(e,e') is identical to (TF(e, e') in 
eq. (32) if ap is replaced by as • From this expression, we 
observe that the corresponding decay length Lc becomes 
very long when the injection current is small and there- 
fore |/l-(x, e)| <^ 1. We thus assume that Lc is longer 
than, or at least of the order of, the length of the su- 
perconducting wire, and adopt the boundary condition 
that /l_ vanishes at each end of the superconducting 
wire. We further assume that the distance between the 
injection junction and each end of the superconductor is 
nearly equal to L, and L ^ d. Under this assumption, 
we approximately obtain 

/L-(a;det,e) = ^7277-l^^^-F^S-(e,^inj) (40) 



Se'^NsAsRiniDs' 



with 



E_(e, Mnj) = Ni{e){fo (e - eMnj) + /o (e + eMnj) 

-2/o(e)). (41) 

li L S> Lc- we must replace L in eq. (40) with L,.- Finally, 
we obtain /F+(2/det) e)- It should be emphasized that Ip 
containing /f+ becomes relevant in the low temperature 
regime where the ordinary contribution /q is exponen- 
tially suppressed due to the opening of the energy gap. 
In this regime, the term with /o(e) — /o(e — eVdct) in 
Pf+ can be neglected. Furthermore, since /t+ can be 
ignored, the dominant contribution to Pf+ arises from 
the term with This indicates that nonequilibrium 
quasiparticles arc created by energy imbalance in the su- 
perconductor. We thus approximate Pf+ as 

n_ / ^ '5(1/ - j/det)-^i(e) . , . 

= 2e^jVFAFfldet (^2) 
Solving eq. (29), we obtain 

/F+(2/det,e) = , hlTJ^'i ^l(^)/L-(^det,6), (43) 
4e iVF^F^det-t>'F 

where the energy relaxation length Ae is given by Ac (e) = 
^DpTcie). Combining this and eq. (40) and noting that 
quasiparticles are populated near the gap edge (i.e., |e| « 
A), we approximately obtain 



/F+(ydet, e) + /F-|-(ydet, -e) 

L 



Ac(A) 



x7Vi(e)S_(e,Mnj) 



3eVde 



(44) 



Prom eqs. (28) and (44), we observe that Jf = at 
Vdet = 0. This reflects the fact that the quasiparti- 
cle distribution /f+ created by the tunneling of energy- 
imbalanced quasiparticles can contribute to the tunnel- 
ing current only when the energy relaxation time for 
/F+(y, e) is different from that for /F+(y, — e).^" That 
is, the energy relaxation process is essential in obtaining 
a nonzero /p. 

We obtain /p and /s by substituting the resulting 
quasiparticle distributions into eqs. (15) and (28). The 



three terms are given as follows: 



X{T) 



If = 



Rdct 

3RsRf LAe(A) 



MVini,T) 



V2(Mnj,T) 



Vdet 
Rdet ' 

A 

ei?det ' 



(45) 
(46) 
(47) 



where = 1(— 1) for the parallel (antiparallel) alignment 
and 



X{T) = 1^ dcNiie) (^-2 



dfoje) 
de 



(48) 



J3{V,T) 



deNi{e){fo{e-eV)+fo{e + eV) 



-2/o(e)), (49) 

1 

MV, T) = ^j deN^ie) (/o (e - eV) - /o (e + eV)) . 

(50) 

The resistances Rs and Rf are defined by Rs = 
XsfPs/Ag and Rp = Asf/Op/^F with the resistivities 
PS = (2e2iVsL's)"^ and pp = {2e^NFDF)-^. It should 
be noted that J3(V,T) and J2iV,T) diverge if eq. (35) 
is adopted as the expression of A^i(e). This unphysical 
divergence does not arise if we adopt a more realistic 
expression of iVi(e), which does not diverges at the gap 
edge. Indeed, the divergence of A^i(e) is actually removed 
if we take account of gap anisotropy, inelastic electron 
scattering or magnetic impurity scattering. We obtain 
Vp and Vap by solving /dct(V^) = for 77 = 1 and 
Idet{Vap) = for 1] = —1, respectively. Substituting the 
resulting expressions and eq. (10) into eq. (1), we finally 
obtain 



R 



spin 



(51) 



with 



Ji{Vini,T) (x(r) 



8-Rinj Rdet 



— 



■MVini,T)y 

(52) 



Note that 7 represents the renormalization of the spin 
signal induced by the transition to the superconducting 
state, and 7 = 1 corresponds to the normal state. 

In the remaining of this section, we briefly consider the 
influence of magnetic impurities. We redefine Tgf as 

14 4 

77 = 3^ + 3^' 

'si '-''SO '->'m 

and introduce the parameter-'^-'^ 



(54) 



which characterizes the relative strength of spin-orbit 
scattering and magnetic impurity scattering. Here, Tgf 
should be regarded as the spin-flip scattering time in the 
normal state. We observe that /3 = — 1 in the absence of 
magnetic impurities and /3 = 1 when spin-orbit scatter- 
ing does not occur. If /? ^ — 1, we must solve eq. (16) to 
obtain the spectral functions. Strictly speaking, eqs. (35) 
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and (36) are not justified in the presence of magnetic 
impurities and the relation iVi (e) — i?2(f) = 1 no longer 
holds exactly. Consequently, /L+(a;det, e) and /L-(a;det, e) 
are modified as 

-Pspina(e)Asf S+(e,Vinj) 



the renormalization factor is reduced to 

^2(/inj,T) 



7 = 



Ji(/inj,T)x(T)- 



(61) 



/L+(a;det,e) 



Se^NsAsRiniDs {Nf{e)-RUe)) 

xe"^^^^, (55) 



L 



S-(e,Mnj) 



/L-(a^det,e) = o„9Ar A o n /Ar2/^ ' p2/^^ ' (^6) 



8e2iVs^sJ?injI)s {Nl{e) - i?i(e)) ' 



Because the T-dependence of J2{Iinj,T)/Ji{Iinj,T) is 
weak, we obtain 7 oc x(^)~'^- This indicates that 7 
behaves as 7 oc c^/'^ . However, because x(T) is expo- 
nentially suppressed with decreasing T, the term with 
J3 (/inj , T) eventually dominates x(^) below a crossover 
temperature Tcross- Below Tcross, we can ignore x{T) in 
eq. (52) and the renormalization factor is reduced to 



where 



a{e) = 



7: 



8i?ini-Rdet 



X2 



^2(/inj,T)J3(/i„j,T) 



(57) 



3RsRf iAe(A) Ji(/i„j,r) 
The crossover temperature is determined by 



(62) 



The parameter a(e) represents the renormalization of the 
spin-flip scattering time on transition to the supercon- 
ducting state. Using eqs. (55) and (56), we can show 
that eqs. (51) and (52) are applicable to this case if J3 
and J2 are replaced by the following expressions. 



X(-^cross) 



3i?S-RF iAe(A) 
SRinjRdet 



J3{V,T) 



NUe) 



foie + eV)-2fo{e)), (58) 



1 



As T is lowered below Tdoss. the injection voltage Vinj 
approaches to A/e. In this situation, the T-dependence 
of J3(/inj, T) becomes weak. Furthermore, we can neglect 
the weak T-dependence of J2(-finj, T)/,7i(/inj, T). Thus, 
we conclude that below T^ross) the renormalization factor 
7 rapidly converges to the value given by 70 = limj'^o 7- 
We can obtain 70 from eq. (62) with T = 0. 

From the above argument, we observe the qualitative 
behavior of iigpin as follows. In the regime of T Tcross; 
the spin signal exponentially increases with decreasing T 



X (./"o (e - eV) - fo (e + eV)). (59) as iispin oc e^/'^ . Below Tcross) however, the spin signal 

^ 7oPspiniise~'^/'^='. We here point 



We here comment on the expression of the spin signal 
presented by Poli et al}^ We note that they ignore the in- 
fluence of nonequilibrium quasiparticles in the ferromag- 
netic metal for detection and therefore the correspond- 
ing term is lacking. This is the significant difference be- 
tween their expression and ours. In addition, they assume 
Vinj "C A/e. Finally, we point out that a(e)-dependence 
is slightly different between them. Indeed, if the factor 
2a+N{E)RN/Ri incq. (4) ofref. 11 is replaced by 20:" 1, 
their expression becomes nearly identical to ours in the 
case of JsiVinjjT) = and Vinj ^ A/e. The reason for 
this difference is not clear. 

4. Discussion 

Let us consider the temperature dependence of the 
renormalization factor 7 under the condition that the 
injection current 7inj is kept constant. We adjust Vinj to 
supply a constant injection current. This means that Vinj 
is determined as a function of T for a given /inj, so we 
rewrite Ji(yinj,T) as Ji(/inj,T) {i = 1,2,3). It should 
be noted here that even though /inj is very small, Vi„j 
approaches to A/e as T ^ 0. 

We focus on the low temperature regime where the T- 
dependence of A can be neglected. In this regime, x(2^) 
behaves as 



X{T) 



2nA 



(60) 



When T is not very low and x{T) is much greater than 
the term with J3(/inj,T) in the denominator of eq. (52), 



converges as Rsj 
out that the behavior of /?spin in the regime of T 3> Tcross 
is qualitatively equivalent to the previous result, eq. (3), 
reported by Takahashi and Maekawa.^ However, our ar- 
gument indicates that the exponential increase of i?spin 
should not be attributed to the increase of spin resistiv- 
ity.^ We simply understand that /?spin increases reflecting 
the suppression of thermally excited quasiparticles in the 
detection junction. 

Let us consider the experimental result reported by 
Poli et al.^^ on the basis of our theoretical framework. 
Particularly, we focus on the convergence of i?spin ob- 
served at low temperatures. They employed the device 
consisting of a superconducting wire of Al and ferromag- 
netic metal wires of Co. Since it has been believed that 
spin-flip scattering in Al is mainly caused by spin-orbit 
interaction, we assume that magnetic impiirity scattering 
is much less relevant than spin-orbit scattering and set 
Tj^^ = 0. We estimate the limiting value 70 of the renor- 
malization factor from eq. (62) with T = and compare 
it with their experimental value. Following rcfs. 10 and 
11, we employ the parameters: /inj = 1 nA, Agf = 1 /zm, 
A = 200 AicV, /?inj = Rdct = 100 m, PS 10 ^ncm, 
As = lOx 150 nm^, Ap = 50x 130 nm^. For the other pa- 
rameters, we assume L = 10 ^m, Dp = 3 x 10~^ m^s"-^, 
pp = 14 /iOcm, Of = 9 X 10^ eV^^. The value of ap 
is estimated by using the relation^^ ap ~ 2/tdTq with 
Td = 385 K and td = 0.4 x 10"^'' s, where To and td are 
the Debye temperature and the phonon scattering time 
at Td, respectively. From these parameters, we obtain 
Ac (A) = 9 iJ,m, Rs = 67 n and Rf = 22 n. The integral 
Ji(/inj,T) does not depend on T and is obtained from 
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eq. (10) as Ji(/i„j,r) = ei?i„j/i„j/A = 0.5. We finally 
consider J2(/inj, T") and Js{Iinj,T) in the limit of T ^ 0. 
The evaluation of these integrals is not simple, so we 
roughly approximate them as J2(iinj,0) = J3(/inj,0) = 
Ji(/inj,T). Substituting these parameters into eq. (62), 
we approximately obtain 70 ~ 10^. This indicates that 
■Rspin below Tcross is a factor of 10^ larger than that in the 
normal state. This is consistent with the experimental re- 
sult which indicates the enhancement of 4 or 5 orders of 
magnitude. Wc estimate the crossover temperature by 
solving eq. (63) with eq. (60) and obtain Tcross ~ 0.1 K. 
This is also consistent with the experimental value of 
Tcross -0.16 K. 

We have shown that nonequilibrium quasiparticles 
with |e| « A are created in the ferromagnetic metal for 
detection by the tunneling of energy-imbalanced quasi- 
particles, and that these quasiparticles contribute to 7det 
in combination with the energy relaxation process due to 
phonon scattering. It should be noted that the energy re- 
laxation of quasiparticles excites phonons near the detec- 
tion junction, leading to the increase of effective temper- 
ature Teff for quasiparticles. If Teg becomes greater than 
Tcross, the convergence of the spin signal is determined 
by this heating effect instead of the convergence mecha- 
nism which we discussed above. The separation of these 
two mechanisms is a future problem for experiments. 

In addition to the heating effect, we have ignored 
charge imbalance. If the injection current is not small, 
we must consider its influences. Charge imbalance pro- 
vides a nearly constant contribution Jq to Is regardless 
of the alignment of magnetizations. Since Iq must be 
cancelled by Iq and Ip to ensure I^et = 0, we expect that 
both Vp and V^p increases with increasing Iq. However, 
if Jq is sufficiently small, the increase of Vp is equiva- 
lent to that of T4p because both /q(Vdet) and IriVdet) 
linearly depends on Vdet when |Vdet| <C A/e. Therefore, 
wc expect that no qualitative change of the spin signal 
appears as long as charge imbalance is not very large. 

In summary, we have studied the transport of spin- 
polarized nonequilibrium quasiparticles in a supercon- 



ducting wire connected by tunnel junctions to two fer- 
romagnetic metal wires, each of which serves as a spin 
injector or detector. We have presented a basic formal- 
ism to determine spin-polarized quasiparticle distribu- 
tions in this system, and obtained an analytical expres- 
sion for the nonlocal spin signal. We have taken ac- 
count of nonequilibrium quasiparticles in the ferromag- 
netic metal for detection, which are created by the tun- 
neling of energy-imbalanced quasiparticles in the super- 
conductor. We have shown that they induce the conver- 
gence of the spin signal at low temperatures. 
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